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Abstract 

In this work I present a complete analysis of proton decay in an 50(10) model previously 
proposed by Dutta, Mimura, and Mohapatra. The 10 , 126 , and 120 Yukawa couplings con¬ 
tributing to fermion masses in this model have well-motivated restrictions on their textures 
intended to give favorable results for proton lifetime as well as a realistic fermion sector without 
the need for fine-tuning and for either type-I or type-II dominance in the neutrino mass matrix. 
I obtain a valid fit for the entire fermion sector for both types of seesaw dominance, including 
013 in good agreement with the most recent data. For the case with type-II seesaw, I find that 
using the Yukawa couplings fixed by the successful fermion sector fit, proton partial lifetime 
limits are satisfied for nearly every pertinent decay mode, even for nearly arbitrary values of the 
triplet Higgs mixing parameters, with only the mode requiring a minor 0(10“^) cancel¬ 

lation in order to satisfy the experimental limit. I also find a maximum lifetime for that mode 
of t{K~^v) ~ 10^® years, which should be tested by eventual experiments. For the type-I seesaw 
case, I find that all six pertinent decay modes of interest are satished for values of the triplet 
mixing parameters giving no major enhancement, with modes other than K~^v easily satisfied 
for arbitrary mixing values, and with a maximum lifetime for K~^v of nearly 10^® years. 
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I. INTRODUCTION 


It has been well-established that a certain class of 50(10) Grand Unihed (GUT) 
models [T] are capable of elegantly solving some of the most prominent problems with 
the Standard Model. One of the more basic yet intriguing features of these models is the 
ability to naturally accommodate a right-handed neutrino, consequently allowing for a 
well-motivated implementation of the seesaw mechanism for neutrino mass m, a long- 
uncontested ansatz that dynamically explains the smallness of (left-handed) neutrino 
masses. The seesaw was originally implemented in the framework of SUSY 50(10) with 
only the 10- and 126-dimensional Higgs multiplets coupling to fermions [2113]; the 126 
vev also plays the role of breaking B — L and triggering the seesaw mechanism, thereby 
creating a deep mathematical connection between the smallness of neutrino masses and 
the other fermion masses. This seemingly-limited yet elegant approach yielded a realistic 
neutrino sector, including an accurate prediction of the value of 6*13 [HI E], long before 
experiments were being done to measure its properties. This so-called “minimal” 50(10) 
model has been explored much more thoroughly over the years by many authors with 
the arrival of precision measurements HHiniE], and it remains a viable predictor of the 
neutrino sector parameters. 

There are however still a number of common difficulties one faces when attempting to 
construct a complete and realistic candidate for unihcation. Furthermore, these difficulties 
are continually being made more severe by new experimental results, which typically 
manifest as new lower bounds on the possible existence of some proposed feature, as part 
of a disheartening streak of null resultsj^ 

Arguably the most problematic feature common to nearly all GUT models arises 
when one examines the lifetime of the proton. In all 50(10) models, heavy 5U(5)-like 
gauge boson exchanges give rise to effective higher-dimensional operators that allow for 

^ The only recent exception to this null trend was the discovery HU of a significantly non-zero (“large”) 
value for the reactor mixing angle ^ 13 . Despite the excitement among experimentalists, the popularity 
at the time of tri-bimaximal mixing models [13] . which prefer 613 ^ 0 , meant the practical elimination 
of nearly an entire class of active 50(10) research. See [HI [13 HU for examples of such models. 
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quark-lepton mixing and, consequently, nonzero probabilities for proton decay 
widths. Furthermore, in SUSY GUT models, several additional decay modes are avail¬ 
able, as each of the GUT-scale Higgs superhelds contains at least one colored Higgs triplet 
that allows for proton decay through exchange of Higgsino superpartners. 

No one yet knows whether protons do in fact decay at all; so far, the lower limit 
on proton lifetime is known to be at least years, and the partial lifetimes for the 

various decay modes have been continually rising through the hndings of experiments im. 
Thus, if any 50(10) model is to be trusted, its prediction for the proton lifetime must be 
at least so high a number. Most minimal SU (5) models have already been virtually ruled 
out by such limits (technically, a few niches of parameter space do still remain). 

There are ways in which the proton lifetime limits can be satished within the frame¬ 
work of a given 50(10) model, but doing so typically requires substantial hne-tuning 
to create rather extreme cancellations (^ 0(10“^)) among the mixing parameters of the 
color-triplet Higgsinos exchanged in the decay. The values of those mixings cannot so 
far be reasonably recognized as more than arbitary free parameters, so to expect multi¬ 
ple instances of very sensitive relationships among them requires putting much faith in 
either unknown dynamics or extremely good luck. Restricting the SUSY vev ratio Vu/vd, 
conventionally parametrized as tan /?, to small values can provide some relief without can¬ 
cellation for Higgsino-mediated decay channels, but such an assumption is still ad hoc and 
may ultimately be inconsistent with experimental hndings; hence it is strongly preferable 
to construct a model which is tractable for any feasible tan/3. 

If however the GUT Yukawas, which are 3x3 matrices in generation space, have 
some key elements naturally small or zero, then extreme cancellations can be largely 
avoided by eliminating most of the dominant contributions to proton decay width. A 
2013 paper by Dutta, Mimura, and Mohapatra [IS] proposed such a Yukawa texture 
for the 50(10) model including a 120 coupling in addition to the 10 and 126 Higgs 
contributions to fermion masses. The authors gave a tentative analysis of mainly leading- 
order relationships between key fermion £t parameters and proton partial lifetimes for a 
neutrino sector arising from a type-H dominant seesaw mechanism; they also provided a 
cursory analysis for a possible type-I solution. 

The work I present in this paper verihes the initial analysis of ref. [IH] by (a) hnding 
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a stable numerical fit to all fermion mass and mixing parameters, including the neutrino 
sector, where values are predictions of the model, and (6) then finding adequately large 
lifetimes for the dominant modes of proton decay using the Yukawa couplings hxed by the 
fermion £t. I grounded the analysis in conservative assumptions, including large tan/3, 
and a comprehensive calculation relying on as few approximations as necessary. The 
pertinent modes of proton decay I checked for sufficient partial lifetimes are p —)■ K~^u, 
K^£~^, 7r~^u, and 7r^£~^, where £ = e,p. I will present solutions for both type-I and II 
seesaw neutrino masses. 

The results not only give satisfactory predictions for the neutrino sector based on 
corresponding charged sector hts, but also adequately predict sufficiently long-lived pro¬ 
tons without relying on the usual large degree of cancellation. Furthermore, I find that 
the ansatz is completely successful in satisfying the proton lifetime limits without any 
need for tuning in the type-I seesaw scenario, while a modest 0(10“^) cancellation is still 
needed in the type-II case to satisfy the partial lifetime limit of the often-problematic 
p —)■ mode. This combination of type-I and II results is precisely contrary to the 

tentative expectations of the authors in [18]; the discrepancy is due mainly to the unex¬ 
pected significance of the effect of rotation to mass basis on the results of the decay width 
calculations, combined with the numerical details of the rotation matrices arising from 
the charged sector mass and CKM fit. 

The paper is organized as follows: in section [n| I give an overview of the S'O(IO) 
superpotential and the fermion mass matrices following from it, followed by the details 
of the Yukawa texture ansatz; in section |III[ I expand further on the model specifics and 
examine general GUT proton-decay logistics in order to derive the needed partial decay 
widths; in section |IV[ I present the fermion sector results of the numerical fitting to the 
measured masses and mixings; in section |V] I present the results of the calculation of the 


important partial lifetimes of the proton; and in section VI, I discuss the implications of 
the results and give my conclusions. 
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II. DETAILS OF THE MODEL 


As mentioned in the introduction, the S'O(IO) model in question has 10-, 126-, and 
120-dimensional Yukawa couplings contributing to fermion masses. The fields are named 
here as H, A, and S, respectively. Thus, the relevant superpotential terms are 

Wy 3 (1) 


where Tj is the 16-dimensional matter spinor containing superhelds of all the SM fermions 
(of one generation) plus the right-handed neutrino, and i is the generation index. 

After the GUT symmetry breaking, SM-type SU{2) doublet representations 
((l,2,—I) -1- c.c) contained in the decompositions of H, A, and S mix with each 
other, and also with the doublets from 126 (which is needed to preserve SUSY invari¬ 
ance) and any additional helds present in the model for GUT-breaking purposes but not 
contributing to fermion masses, such as 210 or 54. These doublets come in pairs with 
conjugate SM quantum numbers, and each Higgs superheld contains one or two pairs. 
The mass matrix A4x> for each set of doublets is determined by the couplings and vevs 
of the GUT-scale superpotential, and so the fields are generally expected to be heavy; 
however, one pair must remain light in order to play the role of the MSSM Higgs doublets 
Hu,d- This need requires the imposition of the fine-tuning condition DetAiv ~ 0 (he., 
^susY ~ 0 when compared to the GUT scale), which can be interpreted as the hxing 
of one parameter in the matrix, conventionally chosen to be the mass of the 10, Mh. 
This choice will have implications for proton decay analysis that I will discuss in the next 
section. In light of this establishment of the MSSM doublets, the effective Dirac fermion 
mass matrices can be written as 


Mu 

Md 

Me 

M., 


h + r2f + r^g 

ri 


tan 13 

ri 


{h + f + g) 

(h - 3f + Ceg) 


tan 13 
h - 3r2f + c^g, 


( 2 ) 


where 1/tan/? takes Vu —t Vd for down-type fields. Each coupling for A = h,f,g is 
related to Xij from eq. ([^ by an absorption of the SUSY vacuum expectation value (vev) 
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Vu and some function of elements of the unitary matrices Ufj, that diagonalize 
These mixings are given in detail with respect to this model in [12] and more generally in 
|21j . but those details are not relevant at this point in the discussion. The coefficients 
and C£ are similarly dehned as functions of those mixings. 

The full neutrino mass matrix is determined by both Majorana mass terms in the 
superpotential and the Dirac mass contribution given in eq. (|^. The light masses can be 
generally given by a combination of the type-I and type-II seesaw mechanisms, involving 
the vevs of both left- and right-handed Majorana terms: 

M, = vl/- (M,„f, (3) 


where are the vevs of the SM-triplet and singlet in 126. The seesaw scale 
(he., RH-neutrino scale) coincides with the B—L breaking scale and is set by u/j; typically 
vl ~ although it is a free parameter of the model in principle. We will separately 

consider cases of type-II {vl term) and type-I {1/vr term) dominance, which can both be 
readily accommodated in this model. Note that the presence of the / coupling in both 
terms intimately connects the neutrino mass matrix properties to those of the charged 
sector matrices, making the model quite predictive. Also note we will consider only normal 
mass hierarchy in this analysis. 

The matrices h and / (with tildes or not) are real and symmetric, and g is pure 
imaginary and anti-symmetric; hence, the Dirac fermion Yukawa couplings are Hermitian 
in general, and their most general forms can be written as 



hii hi2 hi3 


fll fl2 /l3 ^ 

h = 

hl2 h22 ^23 

/ = 

fl2 f22 /23 


hl3 ^23 M y 


y /l3 /23 fsS j 


9 = 


^ 0 gi2 5'i3 ^ 

—gi2 0 g22, 

y ^1713 —92‘i 0 y 


( 4 ) 


M = ^33 ~ mt is singled out to stress its dominance over all other elements. The 
three matrices as written have a total of 15 parameters; taken in combination with vl 
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as well as the vev and mixing ratios r* and q, the model has a total of 21 parameters. 
Correspondingly, there are in principle 22 measurable observables, including all masses, 
mixing angles, and CP violating phases, associated with the physical fermions, although 
the three PMNS phases and one neutrino mass are yet to be observed. Therefore we would 
prefer then to have no more than 18 parameters in the model, and generally speaking fewer 
parameters indicates greater predictability. 

Furthermore, as I will discuss in more detail shortly, the d = 5 effective operators 
that arise in proton decay are ~ XijXki (again A = h,f,g); therefore, increasing the 
number of Xij that are small or zero will increase the number of negligible or vanishing 
contributions to the decay width. This idea was given thorough consideration in [18], and 
the couplings suggested by the authors are as follows: 


h = 


(o \ 

0 

m) 



^ ~ 0 ~ 0 /i3 ^ 

~ 0 /22 /23 ’ 

/l3 /23 fss y 


9 = 


^ 0 gx2 5'13 ^ 

~ 5 ' i 2 0 g2‘i 

— 9 \Z —923 0 y 


(5) 


Note that h is an explicitly rank-1 matrix, with M ~ 0{1)] thus, at first order, the 

10 Higgs contributes to the third generation masses and nothing more. This feature 
has been explored in models demonstrating a discrete flavor symmetry in e.g. [TBl I2U] . 
and may therefore be dynamically motivated. Taking /12 ~ 0 is equivalent to a partial 
diagonalization of /, so it can be done without loss of generality, and the restriction on 
/ii is clearly phenomenologically motivated given the small first-generation masses. As a 
result of these assumptions, the above Yukawa texture should give rise to sufficient proton 
decay lifetimes without the need for the usual extreme cancellations. 

It is further preferred for proton decay that / 13 , 5^12 <C 1, although /13 plays a role 
in setting the size of the reactor neutrino mixing angle ^ 13 , so the above restriction may 
create some tension in the fitting. 

In carrying out the numerical minimization, I will allow /n and /12 to have small 
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but non-vanishing values, 0(10 for the sake of giving accurate first-generation masses 
without creating tension in other elements. The results of that analysis will be discussed 


in section IV, after I discuss the details of calculating proton decay. 


III. DETAILS OF PROTON DECAY 

In addition to the the SM-doublets present in each of the GUT Higgs superfields, 
which contribute to the emergence of Hu,d at the SUSY scale, the GUT fields similarly 
contain SM-type SU{3) color-triplets ((3,1,—|) -f c.c) in their decompositions. These 
fields will also mix after the GUT-scale breaking (again, this mixing includes triplets 
contained in the Higgs helds not contributing to fermion masses). Since there is no 
light triplet analog to found in the low-scale particle spectrum, all of the fields can 
be heavy, although the decoupling of the doublet-triplet behavior is a substantial topic 
itself. The Yukawa potential in eq. ([^ leads to interactions with these helds of the forms 

{QL + U^D^) and hH^ {QQ + which violate baryon or lepton number. The 

helds T'' = are left-handed anti-fermion superhelds. Note e.g. is shorthand 

for the S'f/(2)-doublet contraction . There are similar interaction terms for S 7 - 

and A 7 -; furthermore, two more exotic types of triplets also lead to B- or L-violating 
vertices, (3,1, —-f c.c, which interact with two up-type or two down-type RH singlet 
fermions, and (3, 3, — |) - 1 - c.c, with a pair of LH doublets. 

Exchange of conjugate pairs of any these triplets, through a mass term or interaction 
with a heavy Higgs held such as 54 or 210 , leads to operators that change two quarks 
into a quark and a lepton; this is the numerically dominant mechanism through which a 
proton can decay into a meson and a lepton; corresponding s-channel decays through the 
scalar superpartners of these triplets, as well as s-channel decays through the S'U(5)-like 
gauge bosons X, Y are suppressed by an additional factor of 1/Mu and so are generally 
negligible in comparison]^ Figure [T] shows Feynman diagrams for two examples of the 
operators in question. 


^ The dominant mode in A-boson exchange, p —>■ 7r°e+, may be comparable if the relevant threshold 
corrections are large. 





FIG. 1: Examples of superfield diagrams that lead to proton decay in this model. The hats on the 
couplings indicate mass basis, and the parameters Xi contain the triplet mixing information unique to 
the specific pairing of couplings present in each diagram (see below). 


The Effective Potential 


At energies far below the GUT scale, the triplet helds are integrated ont, giving 
fonr-point effective snperfield operators, which give rise in tnrn to fonr-ferniion operators. 
The corresponding d = 5 snperpotential is 


TVa_b=i — 


^abc 


M- 


T 


/^a./^b^c T I jjC a'pfibjjC c TpC ' 


( 6 ) 


where i,j, fc, Z = 1, 2, 3 are the generation indices and a,b,c= 1,2, 3 are the color indices; 
the SU{2) doublets are contracted pairwise. This potential has AL = 1 in addition to 
AB = 1 and so also has A{B — L) = 0. Mj- ~ Mu is a generic GUT-scale mass for the 
triplets. Note the anti-symmetrization of i, k in the Cr operator; this is the non-vanishing 
contribution in light of the contraction of the color indices. The analogous anti-symmetry 
for the L operator is ambiguous in the current notation, but I will tend to the issue shortly. 

The effective operator coefficients Cijki are of the form 

^ijkl ^O^ij^kl T ^ifijfkl T ^29ij9kl T T ijhkl T ^5fij9kl 


T XQfjij fjf;i -fi -|- x^Qijhki xgfugj}^ -t- xiogngjk 

^ijkl Xghijhj^i Xifijfki X^hij fX4.fijhj^i -\- y5fij9kl T y7hij9kl 

+ y^yikfji + yw9ik9ji- (V 


The couplings h, /, g as written correspond to matter fields in the flavor basis and undergo 
unitary rotations in the change to mass basis, as indicated by the hats on in eq. ([^ 
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above; I will save the details of the change of basis for later in the discnssion. The 
parameters Xi^yi ~ U^,Vjj are elements of the nnitary matrices that diagonalize the 
triplet mass matrix Aij-, or the corresponding matrices for the exotic triplets. Note that 
several identihcations have already been made here: j/o,i = xo,i and j/ 3,4 = —X 3 , 4 ; the 
wonld-be parameters 2 / 2 , 6,8 = 0. Also note that xq ~ Mh ~ 1 is the 10 mass parameter 
hxed by the tnning condition for Mp. The parameters X 9 ,io and 2 / 9,10 correspond to the 
exotic triplets; the indices of those terms are connected in nniqne ways as a resnlt of the 
distinct contractions of helds. 

The left-handed term in eq. (§ can be fnrther expanded by mnltiplying ont the 
donblets as 

WiB., 3 ^ . (8) 

where M is the left-handed nentrino snperheld. Note that the coefficients are sym¬ 
metrized in i, j, as a resnlt of the donblet contractions, and anti-symmetrized in the indices 
of the like-flavor qnarks, again dne to the anti-symmetry of color index contraction, as 
discnssed above for C^. This anti-symmetry will be crncial in restricting the nnmber 
of contribnting channels for decay. Since the symmetrizing of i,j is the same for both 
types of left-handed operators, I will snppress its denotation in fntnre instances to let 
readability favor the less trivial anti-symmetry. 

Dressing the Operators 

As holomorphism of the snperpotential forbids terms like scalar 

boson components of the triplet snperfields, diagrams of the type in Figure can only be 
realized at leading order through conjugate pairs of Higgsino triplet mediators. Thus, in 
component notation, each vertex will be of the form A (p-j-Q Q or similar. Therefore, the 
squarks and sleptons must be “dressed” with gaugino or (SUSY) Higgsino vertices to give 
d = 6 effective operators of the four-fermion form needed for proton decay. Depending on 
the sfermions present, diagrams may in principle be dressed with gluinos, Winos, Binos, 
or Higgsinos. Examples of appropriately-dressed component-field diagrams which give 
proton decay are shown in Figure 

In the following subsections, I will briefly discuss the implications for each type of 
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FIG. 2: Examples of dressed diagrams leading to proton decay in the model. Diagram (a) shows a 
contribution to p —>■ integrating out the triplets gives an effective operator of type C^udue. 

Diagram (b) shows a C^uddiy-type operator contributing to Diagram (c) shows a 

d^-type operator contributing to K^e'l, for Z = 1, 2. Note where more than one field is listed, 
each choice gives a separate contributing channel, except for the dependent exchange of (s O d) in (b). 


dressing and determine which types will contribute leading factors in the proton decay 
width. Note that I will give this discussion in terms of B, W^, and rather than 
and because (a) I will assume a universal mass spectrum for superpartners to satisfy 
FCNC constraints, meaning the mass and flavor eigenstates coincide for the gauge bosons, 
and (6) the mixing of Higgsinos, while not typically negligible, will result in chargino or 
neutralino masses different from Msusy by 0{1) factors, as long as gaugino soft masses 
are relatively small compared to Msusy] since precise values of such masses are insofar 
unknown, and since so many of the SUSY and GUT parameter values needed for the 
decay width calculations are similarly unknown, I will take m^± ~ m^o ~ Msusy = p in 
order to simplify the calculation, especially for computational purposes. 

Gluino Dressing 

Two limitations are readily apparent when considering dressing by gluinos. First, 
the lepton will have to be a fermion leg in the triplet exchange operator, as in Figure [^(b) 
or (c), since a slepton cannot be dressed by a gluino. Second, since SU{3)c interactions 
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are generation-independent, the gluino can only take u ^ u, s ^ s, etc. The latter may 
seem a fairly innocuous idea on its own, but consider that proton decay to a kaon or 
pion will involve operators with one and zero second-generation quarks as external legs, 
respectively, with all others hrst-generation. Taking these two points together with the 
generation-index anti-symmetry of the Cijki operators, which implies that i ^ k for the 
UiDjUkEi operators and j ^ k for the UiDjDkAfi operators, one can see by inspecting a 
dressed diagram that only diagrams with exactly one each oiU,D,S in the triplet operator 
may be successfully dressed by the gluino. This constraint implies that gluino dressing 
can contribute only to p —)■ K~^9 decay mode; furthermore, the absence of UDUE-tjpe 
contributions implies no right-handed channels. 

Taking these constraints into account, and thus looking specihcally at variants of the 
UDSAf operator, there are three independent terms we can write [22], which correspond 
to the dressed diagrams shown in Figure 


Applying the gluino dressing to each term gives us the following sum of four-fermion 
effective operators: 

^abc -F K2{d’’ui){u‘^s'^) -f , (10) 

where the parameters Kq contain factors from the scalar and gluino propagators in the 
loop integral. The scalar propagators are different in general; however, recall that I am 
assuming universal sfermion mass prescription, meaning that all squark masses a equal 
to leading order. In that case, all ks are equal and can be factored out of the brackets. 
The sum left inside the brackets is zero by a Fierz identity for fermion contractions [23] . 
and so the contribution from gluino dressing to the K^9 decay mode vanishes under the 
universal mass assumption. 


Bino Dressing 

As with SU{3)c, U{1)y interactions are also flavor-diagonal; thus, the same con¬ 
straints apply here as in the gluino case, and possible contributions are to the K^9 mode 
^ Each term like is actually the details have been suppressed simply for readability. 
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D 



FIG. 3: Gluino dressings of the d = 5 operator that would contribute to p —>■ 

in the limit of universal squark masses, the three diagrams sum to zero by a Fierz identity. 
NOTE: gluino mass insertions have been omitted from the diagrams for readability. 



FIG. 4: Bino dressings of the d = 5 operator DSN involving a scalar neutrino that would 

contribute to p —>■ K'^Di\ again, in the limit of universal squark masses, the three diagrams sum to zero 
by a Fierz identity. NOTE: Bino mass insertions have been omitted from the diagrams for readability. 


only. 


Looking again at the UDSM operator, for terms in which the neutrino is a fermion 
leg, the argument is analogous to that given for the gluino dressing: the diagrams involved 
are identical to the three in Figure [^except with g ^ B; starting again from expression ([^ 


and applying the Bino dressing, we arrive at an expression similar to (10) but containing 
hypercharge coefficients in addition to the Ka- 


A eabc (A) {K,Y,Ys{u^ui){d’’s^) + K 2 Y^Ys{d%){u‘^s‘^) 


( 11 ) 


however, u,d,sE Qi are all left-handed quarks with Y = 1, so the hypercharge products 
factor out, and again the fermion sum vanishes by the Fierz identity. 

Because leptons carry hypercharge, there are three additional diagrams one should 
include in Figure if dressing instead by the Bino, namely, those involving the scalar 
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neutrino; these diagrams are shown in Figure and the corresponding terms from the 
triplet operator are: 

eabcU^D^S^Ni 3 |(dV)(hV0 + . (12) 


Applying the Bino dressing to each of these terms gives us another sum of four-fermion 
effective operators involving hypercharge: 


B 


+ nn(rf)(sV0}; 


(13) 


this group of terms has a different product of hypercharges from that of ( |IT| , but it 
still has a single common product among the three terms, so we can again factor it out, 
leaving us with yet another vanishing contribution by the Fierz argument. Hence, the 
entire Bino dressing contribution to the mode also vanishes under the universal 

mass assumption. 


IFino Dressing 


As the flavor-diagonal restrictions of the gluino and Bino also apply to the but 
not to the hF^, they must be considered separately. That said, one additional restriction 
applicable in both cases is the ability to interact with only left-handed particles; thus 
there will be no contribution here from the i?-type operators. 

Neutral Wino. As noted, dressing with the is also restricted to UDSM con¬ 
tributions to the mode. The terms to be dressed are the same as those in the Bino 
case, given by expressions ([^ and (12); however, in applying the dressing, we hnd a kink 
in the previous argument: 


fttak (2) {TjTi>(«V,)(<;''!>') + T^THd^UuV) + )(«“/)} 


( 14 ) 


(^) {r=r3(ci‘s'=)(«v,) + nTl(u‘‘s‘){d<-v,) + ri>rJ(uV‘)(sV,)} 


^ ^abc ( ^2 


47r/ 




(15) 
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the negative weak isospins carried by the down-type helds prevent us from using the Fierz 
identity argument. Thus it seems we have hnally found a non-vanishing contribution to 
proton decay, albeit to only this one mode. 

There is something yet to be gained from the Fierz identity in this case: the same 
zero sum we have seen in the previous cases tells us that in each expression here, the 
sum of the two negative terms is equal to the hrst term; furthermore, note that the hnal 


expressions in (14) and (15) are actually identical. Therefore, we can collect the above 
contributions into one expression: 


2 X ^ f(-2){(u'^s‘^)(d%) + 


4 V47r 

^ ^abc 


47r/ 




( 16 ) 


Including the factors from the triplet operator, we can write an operator for the entire 
neutral Wino contribution to K~^P: 


(^) + (fV‘)(sV,)}, (17) 

where the sign cancels with that from the UDDM term in eq. ([^. The details of k will be 
discussed in the next subsection. Note I could have instead written the above expression 
in terms of {(f s^){u°'i'i) alone; I chose the two-operator version because the up-up- and 
down-down-type pairings in the latter option are forbidden in Higgsino and charged Wino 
modes and so are not seen in the calculation otherwise. 

Charged Wino. The assumption of universal mass means that the sfermions are 
simultaneously flavor and mass eigenstates; therefore, the would-be unitary rotation ma¬ 
trix for each is simply the identity, ~ 1. As a result, the unitary matrix present in 
the fermion-sfermion-Wino couplings is not Wkm, but rather the single unitary matrix 
corresponding to the fermion quark rotation. Nonetheless, this rotation allows for the 
mixing of generations at the dressing vertices, and the limitations found on the neutral 
current dressings are not applicable. This is quite crucial since it allows for contributions 
from diagrams with any squark propagator not forbidden by the anti-symmetry of the 
C^ki operator. Proton decay modes involving neutral Kaons or pions, which have uu or 
dd as external quarks, would be intractable without generation mixing. Such mixing will 
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of course come at the expense of suppression from an off-diagonal element in the perti¬ 
nent unitary matrix, which will typically be 0(10“^”^); hence, one can begin to see an 
indication of why the mode is so dominant in the full proton decay width. 

One additional constraint on charged Wino dressing involves the Wino mass inser¬ 
tion. Unlike the gauginos discussed so far, are the antiparticles of each other, rather 
than either being its own antiparticle. As a result, the Wino mass term is of the form 
M^W^W~ in order to involve one and one W~ in the dressing, the two sfermions 
involved must be of opposite SU{2) flavor. As a result, triplet operators of the form udue, 
udue (or the RH equivalents), uddu, and uddu do not contribute. 

Beyond these constraints, the generational freedom of the sfermions leads to nu¬ 
merous contributions to each of the crucial decay modes, K^v, and 

where i = e,fi. In particular the UDUE- and UDDM-tYY>^ operators each contribute 
to each mode through multiple channels. A list of all such contributions would likely be 
overwhelming to the reader no matter how excellent my choices of notation, but one can 
hnd the relevant diagrams in Appendix [A} 

Higgsino Dressing 

When compared to the others, Higgsino dressing is wildly unconstrained. First, the 
low-scale Yukawa couplings governing the fermion-sfermion-Higgsino interactions couple 
a left-handed held to a right-handed one, so clearly the dressing can be applied to both 
C^- and C^-type triplet operators. Also, since charged and neutral Higgsinos couple 
through the same Yukawas, both types of interactions can mix generations, meaning the 
generation-diagonal constraints on the rest of the neutral-current dressings do not apply to 

The only previously-mentioned restriction that does apply is, like the charged Wino, 
the mass term for the SUSY Higgs couples Hu to Hd, so it therefore cannot contribute 
through the triplet operators with sfermions of like SU{2) havor. One remaining minor 
restriction is that we will not see the triplet operator udue dressed by nor uddu dressed 
by because each would result in an outgoing left-handed anti-neutrino. 

One can hnd cases in the literature {e.g. [23]) of Higgsino-dressed contributions being 
counted as negligible when compared to those from the Wino; this is usually because if one 
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exchanges the g 2 Vcabibbo found in a typical dominant Wino contribution for a yfi^, tan f3 
found in a typical dominant Higgsino contribution, the resulting value will be smaller by 
at least a factor of 0(10). Of course one makes several assumptions in such a comparison: 
fi ~ for one, but additionally that (a) tan (3 is small or moderate, and (b) the Cijki 
coefficients are usually of roughly the same magnitude for any combination of 
present. 

For this analysis, though, neither assumption is valid: I have already mentioned that 
I will consider large tan f3 for maximal applicability; furthermore, due to the rank-1 texture 
of the h coupling and the related sparse or hierarchical textures of / and g as shown in 
eq. (j^, many of the Cijki are small or zero, creating large disparities between the values 
from one contribution to the next. This discrepancy from expectation is further enhanced 
by the tendency for the unitary matrices , which give the off-diagonal suppressions at 
the dressing vertices in this model, to individually deviate from the hierarchical structure 

of fo;km- 

To see the extent to which these two properties can lead to surprises in numerical 
dominance, consider that, for example, I hnd 6*^213 ~ ^^ 213 ^ 31 , one might expect that 
U^i ~ Vub and therefore the former term is much larger than the latter, but in fact neither 
assumption is accurate. 

As a result of these model characteristics, I hnd that the dominant contributions 
from Higgsino-dressed diagrams are generally comparable to those from Wino-dressed 
diagrams. This statement further applies to contributions from right-handed operators as 
well. Thus I made no a priori assumptions about which of the C^- or Cfotype Higgsino- 
dressed contributions might be excluded as negligible. 

Because both the operators and the dressing contribute to all of 

the pertinent decay modes, the complete list of channels dressed by the Higgsino is con¬ 
siderably more plentiful than that of the Wino and would again, I suspect, be of no more 
than marginal use to any but the most involved reader. Again though one can hnd all of 
the pertinent diagrams in Appendix 
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Building the Partial Decay Width Formulae 


As I discussed above in the Higgsino dressing subsection, the Yukawa texture seen 
in eq. (|^ leads to (a) unusually extreme variation in the sizes of the Cijki coefficients, de¬ 
pending strongly on the index values present, and (6) textures for the unitary matrices 
which deviate substantially from that of Kkm- The repercussions of these features clearly 
extend beyond affecting the relative size of Wino and Higgsino channel contributions. For 
one, the off-diagonal suppressions f//^, present in most charged Wino diagrams can not 
be dependably approximated as 14^^^A; fortunately, the GUT-scale are hxed by the 
fermion htting, and since the running of such unitary matrices is small, I can simply use 
them at the vertices as reasonable approximations to their low-scale counterparts. 

Another complication due the Yukawa texture is the disturbance of typically useful 
assumptions about which channels dominate the calculation. Such assumptions include 
dominance of Higgsino channels with t, b, f intermediate states or Wino channels oc Va 
or Vcabibbo- In the absence of the validity of any such simplihcation, I am compelled to 
presume that any channel might be a non-negligible contribution to decay width. 

Thus, I initially treated all possible channels as potentially signihcant; however, in 
the interest of saving considerable computational time, I chose an abridged set of contri¬ 
butions to include in my numerical analysis through inspection of tentative calculations, 
although my threshold for inclusion was quite conservative. It seemed to me that conven¬ 
tional methods of keeping only the most dominant terms for calculation might easily lead 
to drastically underestimated decay widths, in that if I exclude ten “negligible” terms 
smaller than leading contributions by a factor of ten, then I have evidently excluded the 
equivalent of a leading contribution. To fully avoid such folly, I used a cutoff of roughly 
1/50 for exclusion, and made cuts on a per-triplet-operator basis, which translates to 
three or four signihcant hgures of precision in the decay widths. 

The Feynman diagrams for all non-vanishing channels of proton decay for the 
K^£~^, TTand modes are catalogued in Appendix [A| 

Calculation of a proton partial decay width can be broken into three distinct parts. 
The hrst part is the evaluation of the “internal”, d = 6 dressed diagrams discussed in 
the previous subsection; each diagram corresponds to an effective operator of the form 
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X qqqi, where X ~ Cijki ... is a numerical coefficient unique to each decay channel. 
Note that here each g is a single quark fermion, not a doublet. The second part is the 
evaluation of a hadronic factor that quantihes the conversion of the three external quarks 
of a dressed diagram-plus one spectator quark-into a proton and a meson. The third 
and final part is the evaluation of the “external” effective diagram for p —)• M£ giving the 
decay width of the proton . I will go through the details of each stage before giving the 
resulting decay width expressions. 

Evaluating the Dressed Operators. The evaluation of one such dressed d = 6 box 
diagram involves calculating the loop integral but no kinematics, because the physical 
particles carrying real momenta here are the proton and the meson, not the quarks. The 
loop factor is not divergent and is of the same general form for every channel; furthermore, 
as the heavy triplets are common to all diagrams and the sfermion masses are assumed 
to be equal, the only factors in the loop that vary from one channel to the next are 
the couplings and masses associated with either the Wino or Higgsino. The remaining 
variation from one diagram to the next depends entirely on the particle flavors, which is 
apparent in the external fermions and encoded in the Cijki coefficients and the unitary 
matrices involved in rotation to mass basis. Thus, I can write the operator for any 
pertinent diagram as a generic Wino- or Higgsino coefficient times one of several flavor- 
specihc sub-operators; the forms of the general operators are 


and 


wher^ 



(18) 


(19) 


^ One might notice that this expression for I{a,b) differs from what is usually given in the literature 
for analogous proton decay expressions; the discrepancy is due to my inclusion of the universal mass 
assumption prior to evaluating the loop integral. 
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and the sub-operators ar^ 

C-' dj) Ui, (dj, C-‘ V,,) 

K = c- e,) <, (7-; (<ij c-‘«) 

'*'=) C'j; t'.l' («^ C-‘ e,,) 

= -^(d C'"’ "') dtili K, (4 C-‘ !■) (20) 

for the (charged) Wino, 

C-^ ei) y]^ (4^ C'“^ 

yji yw 

= -(4 nO 4; 4/ (dr C*-'«") 

<41 = («"^ 4) 4Si]z ylf yw (4 

<^4 = (n"^ ef) a£.,], 4,, 4 (4 n) (21) 

for the charged Higgsino, and 

“4^0 = -(m^ c'“^ 4) r 4'^ (m^^ c-^ 4) 

4^0 = C-^ ei) yll y]! (44 u^ 

Ko = (4 ^mi y'i^ ykl 4) 

4 ) 4 S 1 ]/ 4 Viv {u^ C-^ ei) 

4^ = 4) 4^jk]i yh yjf c-^ 4 ') ( 22 ) 

for the neutral Higgsino, where I have suppressed the color indices everywhere. Again the 
hats on indicate h, /, g are rotated to the mass basis, which I will discuss in detail 
shortly. Note that UDUE and UDDAf operators generally differ by a sign, as do diagrams 
dressed by hu,S latter difference arises from the SU{2) contraction in the 

SUSY Higgs mass term. These sign differences create the potential for natural cancellation 
within the absolute squared sums of interfering diagrams, and even for cancellation of 
entire diagrams with each other in some cases. Also note that the Yukawa couplings are 

® I do not list the neutral Wino operator again here, but looking back at eq. ( [TtI ), we can see that 
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Hermitian in this model, hence the distinction above between and is not relevant 
for this work. 

I utilized two additional observations to simplify the implementation of the above 
operators. First, I took values for the superpartner masses such that /r, <C rriq, which 
imples I{a,b) ~ a/lP‘. Also, because I’m only interested in the combined contribution 
of the three neutrinos, and because the total contribution is the same whether one sums 
over flavor states or mass states, I made the replacement Uii, —)■ 5ui for and took 
1 = 1' ^ 2/fz' = m^i/vd for 

Since the unitary matrices do not appear in the SM (+ neutrino sector) La- 
grangian except in the CKM and PMNS combinations, the non-diagonal SUSY Yukawas 
y^ present in the are not physically determined. Fortunately in our GUT model full 
high-scdXe Yukawas are defined by the completely determined fermion sector. Further¬ 
more, it is known that unitary matrices such as the CKM matrix experience only very 
slight effects due to SUSY renormalization. Thus, since the low-scale masses are of course 
known, I can define good approximations to the SUSY Yukawas needed by using the 
high-scale to rotate the diagonal mass couplings at the proton scale, divided by the 
appropriate vevs: 

(MZ'^)° ul 

where Vu = UwkSin/?, or, in component notation, 

ft" =- (23) 

k 

I can similarly write 

k 

k 

where Vd = v^]^cosf3. Mass values used were taken from the current PDG |23]; light 
masses are run to the 1-GeV scale, top and bottom masses are taken on-shell. Note 
that since the Yukuwa factors always appear in pairs of opposite flavor in the Higgsino 
operators, and since ~ tan (3 for large (3, the Higgsino contributions to proton 

decay ~ for this model. 
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There are generally two distinct mass-basis rotations possible for each of the 
UDUE-, UDDM-, and triplet operators; the difference between the two 

depends on whether the operator is “oriented” (he., in the diagram) such that the lepton 
is a scalar. For a given orientation, a unitary matrix corresponding to the fermionic held 
at one vertex in the triplet operator will rotate every coupling present in pertaining 
to that vertex; an analogous rotation will happen for the other vertex in the operator. For 
example, looking at the channel in Figure j^a), every coupling A^- from present 
at the (f)j- vertex will be rotated by some form of similarly all A).; present at the 
vertex will be rotated by some t/“. The down quark held shown is a mass eigenstate quark 
resulting from unitary the rotation, which we can interpret as a linear combination of ha- 
vor eigenstates: dj = d'^, with j = 1; applying the same thinking to the up quark, we 

can also write = u'^ , with k = 1. To work out the details of the rotations, we can 

start with the d = 5 operator written in terms of havor state^ Xa(ui A“„ d'^)(u'p\'pi q), 
where I have expanded in terms of its component couplings and chosen the indices 
with the malice of forethought; now we can write 



Using the new dehnitions for A, we can see that the rotated coefficient corresponding 
to the expression in eq. ([^ has become 

^ijki Xohijhki -|- Xifijfki x^hij fki 

= x„(;!(7*)y((7;ft)a + x,(f Ul)i0-J)u - xdhU'^U,(U:f)u + ... (24) 

Note that this version of is only valid for U DU E-type operators with this orientation 
in the diagram, namely, those with a scalar e; there is an analogous pair of rotations for 
UDUE with a scalar down and fermionic lepton, as well as two each for UDDM and 
U'^D^U'^E'^, for a total of six possible schemes. 

® Recall the scalars are both mass and flavor eigenstates under the universal mass assumption. Also note 
“A'” is again my name for the second generic coupling, and the prime has nothing to do with basis; I 
will continue to use hats to indicate rotated couplings. 
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From Quarks to Hadrons. As mentioned above, the composite hadrons p and K, n 
(in addition to the lepton) carry physical momenta in the proton decay process, not the 
“external”, “physical” quarks we see in the dressed operators above. Therefore we are in 
need of calculating a factor like (M| {qq)q\p), where M = A', tt is the hnal meson state. 
More explicitly these objects will look like 

{K^\eaUu^s^)Ldl \p) 

\P) 

(7r°| eabciu^d'')Lu'^ji \p) 


Such matrix elements are calculated using either a three-point function (for M, p, and 
the {qq)q operator) on the lattice or chiral Lagrangian methods; in either case, the result 
is determined in part by a scaling parameter (3h defined by (0| {qq)q |p(s)) = (3hPlUp{s), 
where Pl is the left-chiral projection matrix and Up{s) is the Dirac spinor for an incoming 
proton of spin s. In principle fdn is not necessarily the same for cases where the quarks 
have different chiralities, but the values usually differ only in sign, which is irrelevant 
when the entire factor is squared in the decay width expression. 

While lattice methods have advanced significantly since the early years of SUSY 
GUT theory, there is still a substantial amount of uncertainty present in the calculation 
of both jdu and the matrix element factors; some groups have even obtained contradictory 
results when applying the two methods in the same work [2H]. Some more recent works 
{e.g. [25]) using more advanced statistics and larger lattices seem to be converging on 
trustworthy answers, but it is still normal to see results vary by factors of (1/2 - 5) for a 
single decay mode from one method to the next, where the values for the matrix elements 
themselves are 0(10) x (3h- Thus I will simply take the admittedly favorable approach of 
using (M| {qq)q |p(s)) ~ fduPup for all modes. 

It is not uncommon to see values as low as fdn = 0.003 used in other works calculating 
proton decay [2H], but while calculated values have indeed varied as much as (0.003 - 0.65) 
over the years [26], the value is now most commonly found in the range (0.006 - 0.03) 
123 , with a tendency to prefer (3h ~ 0.015, as seen in [26]. Again, I will take a slightly 
optimistic approach and use [3h = 0.008. 
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FIG. 5: Proton decay to a meson and an anti-lepton for; the effective operator vertex contains hadronic 
and renormalization factors as well as the sum of all d = 6 dressed operators contributing to the mode. 


The p —)■ M£ Effective Diagram and the Decay Width of the Proton. Ultimately 
it is a deceptively simple two-body decay that I am calculating, as shown in Figure 
The corresponding decay width can be determined by the usual phase-space integral 
expression: 


1 r d^p 

(27r)3 2UM 


d^p 

(27r)3 2 U, 


(27r)^ S*{pp -pM-Pi)^Yl ^ 

S 


(25) 


where in this case 


^ ^ (|^~|2 +|^-| 2 ) '^\vJ{p^^s)C ( 26 ) 

s s,s' 

The factors and As arise to due the renormalization of the d = 6 dressed operators, 
from Mp to Msusy and Msusy to Mu, respectively; their values have been calculated 
in the literature as Au = 0.4 and Ag = 0.9-1.0 izg. The spinor factor can be evaluated 
with the usual trace methods; in the rest frame of the proton where —Pm = P£ = P, and 
utilizing |p| ^ (which is only marginally valid for the muon but clearly so otherwise), 

the decay width expression simplihes to 


where 


Note that p rNJ 
kaon modes. 


p = Ipl ~ 




1 - 


m 


M 


M 2 


(28) 


Mp/2 for pion modes, but that value is reduced by a factor of ~25% for 
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I now have all the pieces needed to write the working formulae for the partial decay 
widths of the proton. Let me hrst dehne as extended forms of the Cijki by 

= CA (,,)(,«) 

= Cfjqq)(qe) 

= Cf,{qq)(qi), (29) 

SO that these coefficients contain the U-^/2 or factors as well as the Ciju of the 

(30) 


operators in (|20|)-(|22|). Now I can easily translate an operator expression like 

«a;2A 1 
47r / Mj- 


w 




M, 

m 




w 

2 I '-"W ' "W 
g 


into a partial decay width statement, 
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m- 


{AlAsYp\&~ + C 


IV I 2 

•'n' I ’ 


(31) 


without losing either information or readability. Note though there is still a “black-box” 
nature to the C'^ (it was there in the operators as well), in that without specifying 


the generation indices of the external dj^ii quarks, the sums in eqs. (30) and (31) could 
just as easily apply to vr+l'. Furthermore, there are at least several channels present in 
each operator that contribute to any one mode, which are determined uniquely by 
the generations of the internal sfermions in addition to those of the external quarks]^ If 
the reader wishes to examine the decay widths at the full level of detail, he or she should 


utilize these expressions along with the operators in eqs. (20)-(22) and the diagrams in 
Appendix 

All remaining limitations aside, I can now present relatively compact and intelligible 
expressions for the Wino- and Higgsino-dressed partial decay widths of the proton for 
generic mode p —)■ M£: 




M 


w 


m- 


1 (A) 


(.4i21s)= p I 5^ Cj 

A£Me 

EC 


(32) 

(33) 


A&Ai 


^ Indeed I could have defined the coefficients with six indices: thereby creating a means of 

alleviating all degeneracy, but I don’t expect such information-dense objects to be so enlightening to 
readers, especially since for most modes, at least the Higgsino-dressed expression would devolve into 
an entire pageful of terms corresponding to the^g.dividual channels. 



















For the numerical analysis, I used the generic values Mj- = 2 x 10^° GeV, = /x = 
100 GeV, and rriq = 3 TeV. Also, let me repeat here that because of the two SUSY Yukawa 
coupling factors in the Gd, which always come in opposite flavor. 




^wk sin/3 cos/3 


tan^ I3 


■^wk 


Before moving on to the fermion sector £t results, let me remark that because the 
Higgsinos vertices change the chiralities of the outgoing fermions, there can be no inter¬ 


ference between Wino- and Higgsino-dressed diagrams, as suggested by eq. (27); however, 
since diagrams for the right-handed operators have outgoing left-handed fermions by 
the same Higgsino mechanism, diagrams for C^- and G'^-type operators with the same 
external particles of matching chiralities do interfere with each other, and so all such 
contributions to a given mode do in fact go into the same absolute-squared sum factor. 


as suggested by eq. (33). 


IV. FITTING THE FERMION MASS MATRICES 

Diagonalizing the mass matrices given in eq. ([^, with the Yukawa textures shown in 
(|^, gives the GUT-scale fermion masses and mixing angles for a given set of values for the 
mass matrix parameters hij, fij, ri, etc. In order to hnd the best £t to the experimental 
data, I used the Minuit tool library for Python [301 El] to minimize the sum of chi-squares 
for the mass-squared differences Am^i (aka Am^) and Am ‘^2 (^ka and the PMNS 

mixing angles in the neutrino sector as well as the mass eigenvalues and GKM mixing 
angles in the charged-fermion sector. Type-I and type-II seesaw neutrino masses were 
each £t independently, so I report the results for each separately. 

Note that throughout the analysis, I have taken Vu = 117.8 GeV, which is calculated 
with tau/S = 55 and for mn to the GUT scale [22|- The corresponding value for the 
down-type vev is Vd = 2.26 GeV. 

Threshold corrections at the SUSY scale are oc tan/3, and so should be large in 
this analysis [3l|. The most substantial correction is to the bottom quark mass, which is 
dominated by gluino and chargino loop contributions; this correction also induces changes 
to the GKM matrix elements involving the third generation. The explicit forms of these 
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corrections can be seen in a previous work on a related model |T6] . Additionally, smaller 
off-diagonal threshold corrections to the third generation parts oi Aid result in small 
corrections to the down and strange masses as well as further adjustments to the CKM 
elements. All such corrections can be parametrized in the model by 


Ai'd = Md + 


ri 

tan I3 


^ 0 0 

0 0 6Vcb 5 

^ ^^^b J 


( 34 ) 


where Aid is given by eq. Q. If I simply take this augmented form for Aid as part of 
the model input, the 6 parameters are hxed by the mass matrix htting, which results 
in implied constraints on certain SUSY parameters and the mass values that depend on 
them, namely, the Higgs and the light stop and sbottom masses. This entire prescription 
and its implications were considered in detail in [16], and in comparing to that work, 
one can see that for large tan/? and relatively small threshold corrections, the resulting 
constraints on the Higgs and squark masses are less interesting, so I will not consider 
them in more detail for this analysis. 


Fit Results for Type II Seesaw 

If one takes the 126 SM-singlet vev vr ^ lO^^GeV (he., the GUT scale), and the 
triplet vev ul ~ 1 eV, then the type-H contribution {vl term) in eq. ([^ dominates over the 
type-I contribution {vr term) by an average of two orders of magnitude in the neutrino 
mass matrix; therefore eq. ([^ reduces to 

~ vif (35) 

Using this prescription, I hnd a fairly large parameter space for which the sum of chi- 
squares is quite low, although some of the output values, such as 613 and the down and 
bottom masses, are quite sensitive to the variation in the minima. This is problematic 
for 9 i 3 especially, since it is known to high experimental precision [ 12 ] . Tables and 
display the properties of one of the more favorable hts; Table gives the values for the 
adjusted model input parameters, and Table [TT] gives the corresponding output values for 


27 






M (GeV) 

106.6 

ri/ tan/3 

0.014601 

hi (GeV) 

-0.045564 

r2 

0.0090315 

fi2 (GeV) 

0.048871 

n 

1.154 

/i 3 (GeV) 

-0.59148 

Ce 

-2.5342 

/22 (GeV) 

-2.06035 


n/a 

/23 (GeV) 

-1.4013 

6mb (GeV) 

-22.740 

/33 (GeV) 

-1.40644 

<5^6 (GeV) 

1.2237 

gi2 (GeV) 

0.018797 

5Vub (GeV) 

4.2783 

gi3 (GeV) 

-0.92510 



c /23 (GeV) 

-3.8353 




TABLE I; Best fit values for the model parameters at the GUT scale with type-II seesaw. Note that 
Cj/, which appears in the Dirac neutrino mass contribution to the type-I term, is not relevant for type-II. 


the fermion parameters, with experimentally measured values included for comparison. 
Note that the down quark mass is seemingly a bit low, which seems to be a general feature 
in this model, but I will discuss in the next section why this is not a problem. The precise 
value of vl for this £t is 1.316 eV, which I chose to £x the overall neutrino mass scale at 
m 3 ~ 0.05 eV. 

In order to calculate the Cijki proton decay coefhcients, as well as for use in the 
neutrino mass matrix (|^, I needed to determine the “raw” Yukawa couplings, h,f,g, 
from the dimensionful couplings, h,f,g, of the mass matrices given in eq. ([^, which are 
obtained directly from the £t; to do so I need to extract the absorbed vev and doublet 
mixing parameters V]j) mentioned in section]^ There is some freedom in the values 

of those mixing elements from the viewpoint of this predominantly phenomenological 
analysis, but they are constrained by both unitarity and the ratios and q, which have 
been hxed by the fermion fit. Again, see [19] for details, or see |T 6 | for an example of such 












best fit 

exp value 


best fit 

exp value 

mu (MeV) 

0.7172 

' ^-0.15 

Vus 

0.2245 

0.2243 ±0.0016 

me (MeV) 

213.8 


Vub 

0.00326 

0.0032 ± 0.0005 

mt (GeV) 

106.8 


Vcb 

0.0349 

0.0351 ±0.0013 

md (MeV) 

0.8827 


J X 10-5 

2.38 

2.2 ±0.6 

ms (MeV) 

34.04 

29.8ltf 

Amii/Am|2 

0.03065 

0.0309 ±0.0015 

mb (GeV) 

1.209 

1.4210:^8 

013 (°) 

9.057 

8.88 ±0.385 

me (MeV) 

0.3565 

0.3565tO:OOf 

012 (°) 

33.01 

33.5 ±0.8 

m^ (MeV) 

75.297 

75.29l0;?5 

023 (°) 

47.70 

44.1 ±3.06 

mr (GeV) 

1.635 


<5cp (°) 

-7.506 





Ex" 

6.0 



TABLE II; Best fit values for the charged fermion masses, solar-to-atmospheric mass squared ratio, 
and CKM and PMNS mixing parameters for the fit with Type-II seesaw. The Icr experimental values 
are also shown for comparison |31], [H], where masses and mixings are extrapolated to the GUT scale 
using the MSSM renormalization group equations (RGBs). Note that the fit values for the bottom 
quark mass and the GKM mixing parameters involving the third generation shown here include the 

SUSY-threshold corrections 


a calculation. The resulting dimensionless couplings corresponding to this type-II £t are 


/ 


h = 


\ 


1.207 


^ -0.00053748 0.00057649 -0.0069772 ^ 


/ = 


V 


0.00057649 -0.024304 -0.016530 
-0.0069772 -0.016530 -0.0165906 


/ 


/ 


9 = I 


0 0.00033485 -0.016480 

-0.00033485 0 -0.0683214 

0.016480 0.0683214 0 




(36) 


/l2 0, this £t satishes gi 2 , fis <C 1 as is desired for proton 


Note that in addition to fu 
decay. 
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Fit Results for Type I Seesaw 


If one instead takes vr ^ 10^® GeV and vr -C leV, then the type-I contribntion is 
dominant over the type-II contribntion, and eq. ([^ becomes 




( 37 ) 


In this case, initial searches again showed that certain ontpnt parameters were qnite 
sensitive to the inpnt and were often in contention with each other or with the de-facto 
npper bounds on the fij needed for proton decay. In the hrst cluster of minima found 
by the htting, the output values for one or more of charm mass, bottom mass, or 6*23 
was much too small; furthermore, those results came with odd, large tunings of certain 
input parameters, such as Ce,i/ ~ (9(100) or 6 mh > 40 GeV. The addition of a small 
type-II correction to the neutrino matrix led me to a new swath of parameter space, and 
ultimately I found a new cluster of minima that did not require the correction. Table 


III gives the values for the adjusted model input parameters for one such pure type-I £t. 


and Table pV] gives the corresponding output values for the fermion parameters. Fits in 
this swath of parameter space still have Cj, ~ 50 and Sruf, ~ 25 GeV, but this value for 0 ^,, 
while slightly strange, is quite readily accommodated by the doublet mixing parameters, 
and such a value for the largest SUSY threshold correction is actually quite moderate for 


large tan/3. The precise values for the 126 vevs used in this fit are vl = 3.48meV and 
vr = 1.21x10^5 GeV. 

Note also that the top and strange masses are quite a bit lower than in the type-II 
£t; however, note I have also quoted different experimental values with which agreement is 
maintained. The differences here come from an update to the work in |32] in determining 
two-loop MSSM RGBs for fermion masses. The update [33] reports notably lower masses 
for all the quarks at tan/3 = 55 and /i = 2.0 x 10^® GeV, especially for the up, down, 
strange, and top masses, due to updates in initial values and methodology. Hence, one 
should not give the specihc values too much weight in such a £t, and I do not consider 
the reported differences to be significant. This same thinking applies for the type-II down 
mass value in Table m 

Again I need to determine the raw Yukawa couplings for proton decay analysis. The 
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M (GeV) 

76.10 

ri/ tan j3 

0.024701 

/ii (GeV) 

0.010130 

r2 

0.24414 

/12 (GeV) 

-0.089576 

rs 

0.00600 

/i 3 (GeV) 

0.93973 

Ce 

-3.3279 

/22 (GeV) 

0.8659 


45.218 

/23 (GeV) 

1.4884 

5mi) (GeV) 

-28.000 

/33 (GeV) 

3.5495 

dVa (GeV) 

-0.84394 

512 (GeV) 

0.20048 

6Vub (GeV) 

0.51486 

513 (GeV) 

0.05352 



523 (GeV) 

0.35153 




TABLE III: Best fit values for the model parameters at the GUT scale with type-I seesaw. 


resulting couplings corresponding to this type-I fit are 


h = 


\ 


1.6152 


/ = 


/ 


V 


0.0001623 -0.00143525 0.01505699 

-0.00143525 0.01387415 0.02384774 
0.01505699 0.02384774 0.05687217 




g = i 


0 0.0068081 0.0018175 

-0.0068081 0 0.0119376 

-0.0018175 -0.0119376 0 




(38) 


Here, we still see /n ~ 0, but each of / 12 , / 13 , and gi 2 is larger by an order of magnitude 
than in the type-II case, which is thought to be unfavorable for proton decay. At the 
same time, gi^ and 5^23 are smaller by an order of magnitude, so it is not clear that the 
net beneht lost is substantial. In the end, a different distinction will give way to success 
for this type-I £t; I will discuss those details in the next section. 
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best fit 

exp value 


best fit 

exp value 

ruu (MeV) 

0.72155 

‘ ^-0.15 

bus 

0.2240 

0.2243 ±0.0016 

rric (MeV) 

212.2 

210.5l^fi 

Vub 

0.00310 

0.0032 ± 0.0005 

mt (GeV) 

76.97 

80.45l|^* 

Veb 

0.0352 

0.0351 ±0.0013 

md (MeV) 

1.189 

0.930 ± 0.38* 

J X lO"*^ 

2.230 

2.2 ±0.6 

rus (MeV) 

20.81 


Aml^/Aml2 

0.0309 

0.0309 ±0.0015 

mi) (GeV) 

1.278 

1.24 ±0.06* 

0i3 n 

8.828 

8.88 ±0.385 

me (MeV) 

0.3565 

0.3565lO;™f 

0 i 2 n 

33.58 

33.5 ±0.8 

m^ (MeV) 

75.29 

75.29l°;?5 

023 n 

41.76 

44.1 ±3.06 

mr (GeV) 

1.627 


f^CP (°) 

-46.3 





Ex' 

1.75 



TABLE IV: Best fit values for the charged fermion masses, solar-to-atmospheric mass squared ratio, 
and CKM and PMNS mixing parameters for the fit with Type-I seesaw. The Icr experimental values 
are shown [35] (* - [33]), |53|; masses and mixings are extrapolated to the GUT scale using the MSSM 
RGBs. Note that again that pertinent fit values include threshold corrections. 

V. RESULTS OF CALCULATING PROTON PARTIAL LIFETIMES 

In order to give an actual number for any decay width, in addition to choosing 
representative values for the triplet, sfermion, and Wino or Higgsino masses, I also need 
values for the Xi and t/j triplet mixing parameters in order to calculate the Cijki values. 
Recall that the 10 mass parameter xq must be 0{1) to allow the SUSY Higgs helds to 
be light; the remaining mixing parameters are functions of many undetermined GUT- 
scale masses and couplings found in the full superpotential for the heavy Higgs helds, the 
details of which can be seen in |2I]. There are nearly as many of those GUT parameters 
as there are independent xs and |/s, so it is not unreasonable to simply treat the latter as 
free parameters. 

Ideally, one would hnd that the width for any particular mode would be essentially 
independent of those parameter values, he., that for arbitrary choices 0 < |xj|, ||/j| < 1, 
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decay mode 

T exp lower limit (yrs) 

p —)• 

6.0x10^3 

p —)• K^e~^ 

1.0x10^3 

p — )■ Ar°/i+ 

1.3x1033 

p — )• TT~^i> 

2.7x1032 

p — )> 

1.3x103^ 

p — )• 

1.0xl034 


TABLE V: Experimentally determined lower limits [35] on the partial lifetimes of dominant proton 

decay modes considered in this work. 


devoid of unlucky relationships leading to severe enhancements, all mode lifetimes would 
be comfortably clear of the experimentally determined lower limits, given in Table |V} 
The reality is quite bleak in comparison. For a typical GUT model, if the proton decay 
lifetimes can be satished at all, one is required to choose x and y values very carefully such 


that either individual Cs or 


EC-" 


are small through cancellations among terms. These 
tunings may need to be several orders of magnitude in size {e.g., = —C® + 0(10“^)), 

and many such relationships may be needed. 

The Yukawa textures shown in eq. ([^ are intended to naturally suppress the values 
of some crucial C values so that the need for such extreme tuning is alleviated. In order 
to test the ansatz, I “simply” needed to hnd a set of values for the mixing parameters 
yielding partial decay widths that satisfy the experimental constraints; the difficulty in 
determining those values inversely corresponds to success of the ansatz. If the ansatz does 
indeed work optimally, I should be able to choose arbitrary Xi and yi values as suggested 
above. Realistically though, the authors of [18] and I expected some searching for a valid 
region of parameter space to be required. 

To perform that search, I designed a second Python program to hnd maximum 
partial lifetimes based on user-dehned mixing values as well as the raw Yukawa couplings 
hxed by the fermion sector htting. Parameter values are dehned on a per-trial basis for 
any number of trials. I started with the most optimistic case by generating random initial 
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values for Xi and Hi (but Xq ~ 1 fixed), with the decay width for minimized by 

adjusting those values in each trial. The minimization was again performed using the 
Minuit tool library. 

The search based on fully random initial values was unsuccessful, in that the 
mode lifetime consistently fell in the year-range for the type-II solution and was 

typically ~ 1x10^^ years for the type-I caseat the same time however all five other modes 
in question were usually near or above their respective limits for those same arbitrary 
mixing values. Hence it was clear that, even with the mode failure, that the ansatz 
was having the desired effect to some extent. Also, note that this type-I solution for 

was short of the limit by only about a factor of five. This is surprising since the 


type-I-based Yukawas reported in eq. (38) fell short of meeting the ansatz criteria. Given 


the differing behaviors of the two solutions, I will report the remaining details in separate 
subsections once again. 

Proton Partial Lifetimes for Type II Seesaw 


To further explore the properties of the “default behavior” of the lifetime values 
in the model, I considered the case in which xq ~ 1 and all other Xj and Ui are set to 
zero; one can see this case as defining a baseline for the partial lifetimes, in that any xq 
terms in the Gs not suppressed by the Yukawa textures are necessarily large, and whereas 
problematic contributions from some other Xk with k ^ 0 may be suppressed simply by 
setting Xfc 1, the xq contributions can be mitigated only through cancellation. 

The corresponding baseline lifetimes for the dominant modes in the type-II case 


are given in Table VI One can see that the K~^P mode decay width must be lowered 
by two orders of magnitude through cancellation of Xq terms by the others. Since it is 
|C|^ that appears in the decay width expressions, the needed cancellation amounts to an 
0(10“^) tuning among the C"^ factors. Furthermore, as it would be equally unnatural 
to see Xfc -C 1 for all /c 7 ^ 0, one should expect 0{1) cancellations to be present anyway; 
therefore, the needed “tuning” is little more than a very ordinary restriction of parameter 
® The Minuit tool used, Migrad, works using a local gradient-based algorithm, so that in large parameter 
spaces, initial values are crucial in locating global minima. 
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decay mode 

baseline for r (yrs) 

baseline in ref. [H] (yrs) 

p —7- 

8.29 xlO^^ 

6.38x10^3 

p —)• K^e'^ 

9.73x1034 

2.52x 103° 

p —^ K^pL^ 

5.68x1033 

6.15x1029 

P —)• TT^U 

4.25x1033 

4.45x1029 

p —)• 

1.08x1036 

3.90x103° 

p —>■ 

6.45x1034 

6.00x1029 


TABLE VI: Hypothetical baseline partial lifetimes determined using type-II solution Yukawas and 
xo = 0.95 with all other Xi,yi = 0. For comparison, I give the analogous results for calculation using 
type-II Yukawas from the 2010 paper by Alterelli and Blankenburg m, which use general Yukawa 
texture. Note in comparing with Table [v| that for our model, only the mode fails to satisfy the 
lower limit, while all modes are well below the limits for the model in El- 


space. 

In order to elucidate the significance of the improvement created by the Yukawa 
ansatz, let us consider the outcome of this baseline calculation for a case with more general 
Yukawa texture. The model from a 2010 paper by G. Altarelli and G. Blankenburg [Hj 
has the same 10 - 126-120 Yukawa structure but with general h and g as in eq. (|^ and a 
tri-bimaximal / having no hierarchical texture]^ Using the parameters reported to give 
a successful fermion £t in the work (see footnote), I obtain the baseline results shown in 
the hnal column of Table lYl One can see here that lifetimes for all modes are far below 
the experimental limits, by factors of 0(10^”®); hence cancellation among the G'^ factors 
must be 0(10“^”^). Such sensitive relationships among these factors are considerably more 
restrictive than the result from our model, and, in the absence of some new symmetry, 
there is no good explanation for those restrictions. 

In order to locate an area of mixing parameter space which yields a sufficient 
lifetime, I wrote a supplementary Mathematica code to search for minima among strongly 

® This model has already been ruled out due to 613 ^ 6-7° typical of tri-bimaximal models. 
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r(p ^ K+u) (yr) 


t(p K+v) (yr) 



FIG. 6: Comparisons of K^v partial lifetime to those of other dominant modes in the model, and that 
lifetime as a function of the 10 mass parameter xq, for the type-II case. Note the unsurprising 

preference for smaller xq. 


abridged versions of | C~ + C~ | and | | that contribnte to the decay widthj^ Specih- 

cally I started with xq terms only, corresponding to the baseline case, and then iteratively 
added back the largest contribntions one by one while readjnsting the initial valnes each 
time. Once all of the most important terms were present, I took the resnlting mixing 
parameters as my initial valnes in the Python code. The resnlting minimization gave a 
large percentage of trials with all six modes exceeding the lifetime bonnds. 

With an allowed region of parameter space fonnd, I expanded my searches to inclnde 
a wider range of valnes for xq. Using six different “seeds” for parameter valnes, all of which 
give every mode snfficient with t[K^v) ronghly twice the experimental bonnd, I created 

and CF[ cancel identically for all contributing channels of both the K^v and modes. 
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r(p —^ A'®e+) (yr) t{p —>• vr+i^) (yr) 


FIG. 7: Comparisons of partial lifetimes among highly-correlated sub-dominant modes in the model 

for the type-II case. 

a large number of trials for which the initial values were distributed normally around the 
seed values and with large standard deviations. The resulting data for such a search is 
shown in scatter plots below. Figure gives the relationships between the iC+z/ mode 
and other representative modes and also the distribution of lifetime for varying 

Xf). Figure shows the relationships between other more closely correlated modes for 
completeness. 

Note the strong correlation between vr+P and which are related by isospin, and 

the extreme correlation between and The latter is due to a manifestation of 

the hierarchical nature of the Yukawas in the Cijki-, as well as minor features such /n / 12 ; 
similar structure is present in the and U-^, which tend to also have 11 ~ 12 or 11 ^ 12; 
these properties result in a straightforward scaling under the replacement / : 1 —)■ 2. 
Furthermore, the same relationship is present between vr^e’'' and . These relationships 
imply that the remaining plots I omitted differ only trivially from the representatives 
present. 

I also performed simple scans in search of a maximum value for r(iF+P), as well as 
taking note of any especially large values in the previous searches. While there does not 
seem to be any analytically-enforced maximum present in the model, I did consistently hnd 
that r > 10^® years was extremely rare, and I never saw a value higher than ~ 6x 10^® yr. 
Given those hndings, combined with the apparent smallness of the swath of parameter 
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space yielding the above results and the low likelihood of a more global minimum based 
on my search methods, I believe that ^ 10^® yr is statistically infeasible in this 

model for type-II seesaw. If such a value does exist, it is likely contained in a vanishingly 
small area of allowed parameter space and accomplished through truly extreme tuning. 
Therefore I will take 10^® years as a de facto upper limit on t{K^V) for the type-II case, 
which will not be accessible by Hyper-K and similar experiments [361 EZ] in the near 
future, but should nonetheless allow the model to be tested eventually. 

The other modes of course have similar limits, but it would seem that all the others 
are substantially higher and thus either far beyond the reach of the forthcoming exper¬ 
iments or beyond the contributions from gauge boson exchange, if not both, with the 
possible exception of r(7r■*■//), which is rather highly correlated with in this model. 
Determining that value is tricky though because if I simply maximize the tt+z/ mode, then 
the mode will be below its bound; thus, there is some question as to how one dehnes 
the maximization. 


Proton Partial Lifetimes for Type I Seesaw 


I begin again by examining the same baseline case for the partial lifetimes, with 
Xq ~ 1 and all other x,, Ui = 0. The resulting values for the dominant modes in the type-I 


case are given in Table VII Here we see a much more favorable situation, in that even the 
mode decay width is sufficient, and in fact the other modes exceed the bounds by 
2-4 orders of magnitude. Hence we expect that virtually all solutions will be adequate for 
modes other than K~^9, and as long as there is no enhancement due to (de)tuning among 
the factors, that mode will be adequate as well. 

This is of course a remarkable improvement over traditional models, yet it seems 
to contradict our expectations given then properties of the fermion £t. Why then is the 
model successful? There are two primary reasons, both of which are quite subtle. The 


first reason is that the smaller values for gi^ and 5^23 seen in eq. (38) do in fact improve 


the situation, as I suggested, while the larger /12 and gi 2 seem to have less impact. Since 
M {h^^) is such an extremely dominant factor in the Yukawas, it is generally the case that 
contributions involving third generation are larger and more important than the others. 
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decay mode 

baseline for r (yrs) 

p —7- 

7.87x10^5 

p —)■ K^e~^ 

5.93x10^5 

p —)• 

2.45x1035 

p — >■ 'K~^V 

2.37x1036 

p — >■ 

6.11x1038 

p — )• 7r°/i+ 

2.27x1038 


TABLE VII: Hypothetical baseline partial lifetimes determined using type-1 solution Yukawas and 
xo = 0.95 with all other Xi, yi = 0. Note in comparing with Table [V| that all modes satisfy the lower 
limits, and most do so by several orders of magnitude. 


The second reason is even more nnexpected, to the point that it was not even 
examined in the preceding works on this ansatz. The nnitary matrices for the charged 
fermions are generally ~ 1, just as one would expect, given the texture of CKM. This 
model is no exception, with off-diagonal terms generally 0(10“^”^); however, with such 
sparse or hierarchical (flavor basis) Yukawas due to the ansatz, these “small” off-diagonal 
elements lead to “small” rotations of h, /, g resulting in relatively substantial chanes to 
the textures of h, /, g. Especially noteworthy are the changes in h —)• h, where some 
previously-zero off-diagonal elements are replaced by the same 0(10“^”^) values seen in 
the UL 

In light of the surprising non-triviality of the basis rotations, if we compare for 
the type-I case: 


Uu 


^ 0.994 -0.1085+ 0.0057f 0.00298 + ^ 

0.1084+ 0.0057Z 0.994 0.0047 + lO-'^z 

^-0.0035 - 10-5* -0.0044 + 10-5i 0.99998 ^ 


= 


0.967 -0.1087 + 0.2309* 0.00175 + 0.001175* 

0.1086 + 0.2308* 0.966 0.03935 + 0.00690* 

-0.0076 - 0.0072i -0.0381 + 0.00613* 0.9992 


A 


(39) 
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to those for the type-II case: 


/ 


= 


0.972 0.2098 - 0.1044i - O.OlOz 

-0.210 - 0.1043* 0.971 -0.00012 - 0.0414i 

-0.0043 -0.001* -0.001 - 0.0423* 0.999 




= 


0.9998 0.00633 - 0.0095* 0.00765 - 0.01117* 

-0.00708 - 0.0095* 0.9983 0.03386 - 0.04514i 

-0.00785 - 0.01054i -0.03401 - 0.04514i 0.9983 




(40) 


we see that the off-diagonal entries are the same size or smaller for the type-I case in every 
entry except 11 ^ 2 ', furthermore, several of the elements involving the third generation 
are smaller by an order of magnitude. These differences may seem rather benign, but 
in fact each of these slightly suppressed values individually translates into a factor of 
10 suppression in most of the dominant Cs, which all tend to involve third generation 
elements. In some cases two or even three such suppressions may affect a single C"^ factor. 
The squaring of factors in the decay width then gives suppressions of generally 2-4 orders 
of magnitude in the lifetimes, which is precisely what one can see when comparing Tables 

lyTiandivrn 

Due to the more favorable circumstances, I was able to locate an allowed region of 
parameter space for type-I simply by running a large number of trials with the type-II 
parameter seeds. I repeated the process of expanding the range of Xq by again choosing 
hve seeds that gave every mode sufficient and roughly twice the experimental 

bound, and I again used those seeds to create scatter plots for a large number of trials. 
Figure gives the relationships between the mode and other representative modes 
and the distribution of t{K^V) as a function of xq, and Figure shows the relationships 
between other more closely related modes. Note the bifurcation of the solution set in each 
plot; I have not yet been able to discover the cause of this behavior. 

Again I performed scans to determine a statistical upper bound for the value of 
t{K^u) in the model. I consistently found that r > 10^^ years was rare and did not see 
a value higher than ~ 3 x 10^^ yr. Given those hndings, I suspect that the de facto upper 
limit on for the type-II case is slightly lower than 10^® years for the type-I seesaw 
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t{p —>■ K^v) (yr) t(p —>■ K^v) (yr) 



t[p ^ K^v) (yr) ° 

FIG. 8: Comparisons of K^v partial lifetime to those of other dominant modes in the model, and that 
lifetime as a function of the 10 mass parameter xq, for the type-I case. Note the unsurprising preference 

for smaller Xq- 

case. Such a value is clearly out of reach of Hyper-K and other imminent experiments. 
Note that as values for the neutral Kaon and pion lifetimes often exceeded 10^® years in 
my hndings involving iC+z/ minimization, the upper limits for those modes are surely 
sub-dominant to gauge exchange as well as out of reach of experiments and so not of 
interest. 

VI. CONCLUSION 

In this work I have presented a full analysis of the nature of proton decay in an 
50(10) model that has 10 , 126 , and 120 Yukawa couplings with restricted textures 
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r(p —>• _ftr®e+) (yr) t(p —>• tt^v) (yr) 

FIG. 9: Comparisons of partial lifetimes among highly-correlated sub-dominant modes in the model 

for the type-I case. 

intended to naturally give favorable results for proton lifetime as well as a realistic fermion 
sector. The model is capable of supporting either type-I or type-II dominance in the 
neutrino mass matrix, and I have analyzed both types throughout. Using, numerical 
minimization of chi-squares, I was able to obtain successful hts for all fermion sector 
parameters, including the 6*13 reactor mixing angle, and for both seesaw types. Using the 
Yukawa couplings hxed by those fermion sector hts as input, I then searched the parameter 
space of the heavy triplet Higgs sector mixing for areas yielding adequate partial lifetimes, 
again using numerical minimization to optimize results. For the case with type-II seesaw, 
I found that lifetime limits for hve of the six decay modes of interest are satished for 
nearly arbitrary values of the triplet mixing parameters, with an especially mild 0 ( 10 “^) 
cancellation required in order to satisfy the limit for the mode. Additionally, I 

deduced that partial lifetime values of t{K^v) ^ 10 ^® years are vanishingly unlikely in 
the model, implying the value can be taken as a de facto lifetime for the mode, which makes 
the model ultimately testable. For the case with type-I seesaw, I found that limits for all 
six decay modes of interest are satished for values of the triplet mixing parameters that do 
not result in substantial enhancement, with limits for modes other than satished for 
nearly arbitrary parameter values; furthermore, I deduced a statistical maximum lifetime 
for of just under 10^® years. Given these results, I conclude that the well-motivated 
Yukawa texture ansatz proposed by Dutta, Mimura, and Mohapatra is a phenomenological 
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success, capable of suppressing proton decay without the usual need for cancellation and 
without compromising any aspect of the corresponding fermion mass spectrum. 
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Appendix A: 

Feynman Diagrams for All d = 6 Operators Contributing to Proton Decay 

Channels for p — y i,l = 1,2,3; (pj- is the Higgsino component of a heavy color- 

triplet Higgs superfield, ^ = H, A, S 
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Channels for p —>■ 7T^£~^ . j = 1,2,3; ^ = 1,2 (o £ = or for diagrams including /' 


instead ? = 1, 2, 3 and /' = 1, 2 
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Channels for p — y i,l = 1,2,3; parentheses indicate coupled choices; the absence 

of a diagrams containing udue dressed by and uddO dressed by is due to resulting external 
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Channels fov p —y . j = 1,2,3; / = 1,2 (o £ = e,ii), or for diagrams including 

instead ? = 1, 2, 3 and /' = 1, 2 
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